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' Abstract 

An imperfect cosmic fluid with energy flux is analyzed. Even though 
^ | its energy density p is positive, the pressure p = — p due to the fact 

6J[) that the metric is asymptotically de Sitter. The kinematical quantities 

for a nongeodesic congruence are computed. The scalar expansion is time 
independent but divergent at the singularity r = 2m. Far from the central 
mass m and for a cosmic time t « H~ , the heat flux q does not depend 
on Newton's constant G. 
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1 Introduction 

The exterior Schwarzschild solution of Einstein's field equations describes the 
stationary gravitational field around a spherically symmetric star or black hole 
(BH). Such asymptotically flat solutions apply to a region outside the central 
object, far from other compact sources. 

A 2nd set of solutions are the cosmological ones that model the behaviour 
of the universe at its largest scale pQ as, for example, the homogeneous and 
isotropic Friedmann - Robertson - Walker (FRW) metric on which the standard 
model is based. The purpose is to combine both classes of solutions and to find 
exact solutions for the gravitational field of a compact object embedded in a 
cosmological background Q]. One means the solutions describe BHs which are 
asymptotically flat and, in addition, dynamical [2J |31 IH 03 [3] . A complication 
arises when a BH is "immersed" in a cosmological background (excepting the 
static de Sitter case): the phenomenon of accreation will be present. Namely, 
a radial flow of energy of the cosmic fluid tends to infall forward the central 
object [TJ[2j[ai8]. 

McVittie [9], in his well known paper, has taken into account the effect of 
the cosmic expansion on local systems. To forbid the accretion of the cosmic 
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fluid into the central object he imposed a special relation between the (time 
dependent) mass of the object and the scale factor. Since in this case the 
stress tensor has no non-diagonal components, there is no radiative heat flow. 
However, the physical meaning of the McVittie spacetime is not completely clear 
and is still under debate. 

Another exact solution that models a spherical object in a cosmological 
context was proposed by Sultana and Dyer |10j . Their metric is conformally 
equivalent to the exterior Schwarzschild metric and the cosmic matter is com- 
posed of two non - interacting perfect fluids, one being pressureless dust and the 
other - a null fluid. In addition, the Sultana - Dyer spacetime docs not belong 
to the class of McVittie models (it is not spatially - Ricci isotropic pQ). 
Throughout the paper we use the geometrical units G = c = 1. 



2 Time dependent conformal Schwarzschild met- 
ric 



Let us begin with a general geometry 

2m, ,o 2m 



ds z = a\t) 
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(2.1) 



where a(t) is a positive function of the time t, m is the BH mass (or the mass 
of a spherically symmetric compact object) and r > 2m. When m = 0, (2.1) 
gives us a conformally flat spacetime. In contrast, when a(t) = 1 we obtain the 
exterior Schwarzschild metric around the isolated mass m. 

With the help of the GrTensor II package and software package Maple, the 
following expressions for the components of the Einstein's tensor G\ are ob- 
tained : 
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where a = da/dt, the Latin indices run from to 3 and the coordinates are 
in order (t,r,9,4>). We observe that in (2.2) there are nonzero non - diagonal 
components. That means a flux of radial energy will be present in the system. 
Therefore, the cosmic fluid must be imperfect [3] and the corresponding energy 
- momentum tensor given by the Einstein equations G^ = SnT^ should have 
terms dependent on the heat flux vector. It is worth noting that the scalar 
curvature 

6a 

a 3 ( — 2in) 

is singular at r = 2m, a property valid also for the McVittie spacetime [3], 
where the mass m is time dependent. In addition, the metric (2.1) has another 



= TT, 2mT (2-3) 
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(expected) singularity at r = 0, where the Kretschmann scalar R ahcd R a b c d di- 
verges. As far as the Weyl tensor is concerned, all its components are of the form 
a 2 (t) C^ d w , where C^ d w corresponds to the Schwarzschild metric. The Weyl 
tensor is, of course, vanishing when m = because the metric is conformally 
flat in that case. 



3 Imperfect fluid stress tensor 

Let us choose now the following congruence of observers 

1 



,0,0,0 | (3.1) 

with u a the velocity vector field associated to "static" observer (dr/dr = 0, u a u a - 
— 1, t being here the proper time). We assume the stress tensor is given by the 
following expression 

Tab = (p + P)u a Ub + P9ab + u a q b + u b q ai (3.2) 

where u a q a — 0. The energy density p, the pressure p and the heat flux q a 
acquire the form 

p = T ab u a u\ p = T r r =T$ = T*, q a = -T b a u b -pu a . (3.3) 
Using now Eqs. (2.2), one obtains 

3d 2 
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(3.4) 



Let us consider the case of an expanding universe (d > 0). Hence, q 1 = q r < 0, 
namely a radial inflow will be present. In addition, the parameters p,p and 
q r are divergent at r = 2m, a consequence of the singularity of the curvature 
invariant scalars there. 

It is worth to compute the total radial energy flowing accross a r — const. 
hypersurface E [11] for some time interval At = ti — t\ 

,b 




W = J T\u b n a ^J-h dt d6 d4> (3.5) 

where n a = (0, ^Jl — ^p/a, 0, 0) is the unit spacelike normal vector to S, h ab — 

9ab — n a n b is the induced metric on S and h = det(h a b)- By means of T = 
md/4-7ra 3 r 2 from (2.2) and u a from (3.1), we get 

to f* 1 , s mAa(t) , . 

W = . / a{t)dt = V ' (3.6) 
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where Aa(t) — afa) — a(t\). For an observer located at r >> 2m, we have 
W = mAa(t), i. e. the gain in the physical mass of the compact object. In other 
words, all the accreted energy contributes to the increase of the object mass. 
We of course have neglected here the back reaction effect, keeping m = const, in 
the metric (2.1). Moreover, W grows indefinitely when r — > 2m. If At = <2 — ti 
is small, so is Aa and, therefore, W is also small excepting near the horizon 
r = 2m. 

The congruence (3.1) is not geodesic as could be seen from the expression of 
the acceleration vector 



a b = (0 



(0, J£_, 0, 0) , (3.7) 



with y/aFob = m/ar 2 yj\ — As far as the other kinematical quantities is 
concerned, the shear and the vorticity tensors of the congruence are vanishing 
(due to the spherical symmetry) while the expansion scalar is given by 

9 = V a u a = 3a (3.8) 



If we formally compute the surface gravity on the horizon r = 2m, one obtains 

i—r — . ml 

{2m) z 4m 

as for a static BH. It is, however, problematic whether one might define a surface 
gravity in our nonstatic spacetime and, in addition, at r = 2m the metric is 
singular. 



4 Time dependent Schwarzschild - de Sitter ge- 
ometry 

Our next task is to take into consideration a particular expression of a(t) and 
analyze that case in detail. A situation not studied till now, to our knowledge, is 
the Schwarzschild BH embedded in a de Sitter universe, written in conformally 
fiat coordinates. We prefer a spacetime conformal to the Schwarzschild one in 
order to preserve the null trajectories and the causal horizon. 

Let us firstly consider the de Sitter line element written in comoving coor- 
dinates 

ds 2 = ~dP + e 2Hi {dr 2 + r 2 dfl 2 ), (4.1) 

where H = l/Ris the Hubble constant, R is the radius of the de Sitter universe 
(from the static form of the metric, related to the cosmological constant A, with 
1/R 2 = A/3, i is the comoving time and dil 2 = d6 2 + sin 2 9d<fr 2 represents the 
metric on the unit two-sphere. By means of the transformation [TJ] 

x ] = ~e- m , (t>0, -l<r?<0), (4.2) 
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the metric (4.1) may be put in a conformally flat form, 

ds 2 = -^(-drj 2 + dr 2 +r 2 dn 2 ), (4.3) 

r\ being the conformal time. We intend to study a BH embedded in the above 
spacetime, namely in the geometry 



ds 2 ' 



, 2m. , 2 dr 2 2 
.(!_ ) dr f+ z-+r 2 dQ 2 



(4.4) 



for to find the influence of an expanding universe on local physics. We have now 
O'iv) — l/^M which leads to 

d a 1 a 1 2 .. . 2 3 

fl= ^ = W a=N' a = "W >0 ' 2aa_a = W (4 ' 5) 

We notice that d > 0, that is the spacetime is accelerating, as it is confirmed 
for our universe. For the above particular value for 0(77) we have from (2.2) for 
the components of the stress tensor 

2mH 2 r) _ SttT^ 
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r 2(J _ 2m\2 M _ 2m\2 

(4.6) 



87T/9 = -8ttT = -8irp = 



2 2m 



We see that the equation of state of the fluid is p = —p, irrespective of the value 
of m. Moreover, they have no a time dependence. For r >> 2m (far from the 
central mass), p acquires the well-known expression 3H 2 /8w for a spatially flat 
FRW universe. Even though the energy density of the fluid is positive and the 
weak energy condition holds, the strong energy condition is not satisfied. We 
have, indeed 

1 „„\ a b 



Tab - ^9abT^j u a u b = -p < 0, (4.7) 

a basic property of the dark energy (DE), due to the negative pressure. 
The kinematical quantities a b , 6 and q a become now 

h / mH 2 n 2 \ f—r — mHln 
a b =[0, _L, 0, 0), Vo^= — r= 



\ _ 2m 



9 = 3H , ^=[o, ,0,0 



(4.8) 



In addition, (4.8) yields 



47rr 2 A /l- ^ 



6 s «-v e = o, q = yftT a = (4-9) 
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We see that, in spite of the time dependance of the metric (4.4), is a function 
on the radial coordinate only. One means that for an observer sitting at r — 
const, the congruence has a constant expansion. It becomes, however, divergent 
when r — ¥ 2m because of the singularity. 

An interesting property of the spatial vector field q a that represents the 
current density of heat is its independence of the Newton constant G. Far from 
the constant mass and for the cosmic time i << H^ 1 (or from (4.2), H\n\ ~ 1), 
the expression of the heat flow becomes (we introduce here the fundamental 
constants G and c) 

c 4 Gm H mc 2 H 
Q = G = l^ 5 "' ( ' 

Let us find, as an example, the ingoing heat flux on the Sun surface, using 
the above approximation. With Ms ~ 2xl0 33 <?, Rs ~ 7x10 10 cto and H m 
2.3xl0~ 18 s~ 1 , we get qs ~ 6.2xl0 13 ergs/ cm 2 s — 6.88xl0 -8 g/ cm 2 s, or a power 
of 1.4xl0 44 ergs/year. 



5 Misner - Sharp mass 

For to interpret the spacetime (4.4) as a model for an inhomogeneity in a FRW 
universe (spatially flat in our situation) it is useful to compute the quasi-local 
Misner - Sharp (MS) mass M [13[ I14| I15[ 116) . with its Weyl and Ricci parts, 
namely Mw and Mr, respectively. It will help to detect localized sources of 
gravity . The mass M is obtained from 

where R. a = dR/dx a is the areal (physical) radius which in our case is given by 
R{t,r) = ra(t) (see 2.1). Eq. (5.1) yields 

■ 2 3 

M (t, r) = ma H ° T „ - (5.2) 

We identify the 1st term from the r.h.s. with the Weyl mass Mw and the 2nd 
term with the Ricci mass Mr that gives the contribution of the fluid to the MS 
energy. We note that Mw is the gravitational energy and ma(t) represents the 
physical mass. 

Let us now apply the formula (5.2) for the spacetime (4.4). One obtains 

m v 

M ^=m + mm^W) (5 - 3) 

The fact that Mw = vn/H\r}\ grows from m to infinity when n varies from —1/H 
to zero is interpreted as saying that energy is accreated from the ambient matter 
onto the central object, as one sees from the expression of q. Ew comes from 
the Weyl curvature and may be identified with the gravitational mass of the 
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central object. The Ricci part of the Misner - Sharp mass may be expressed in 
terms of the density p from (4.6) and the areal radius, to give Mr = p 4irR 3 /3, 
as expected (see, for instance, pQ). 

We wish now to evaluate the MS mass at t = or \rj\ = l/H, for r >> 2m. 
While Mw ~ m, the Ricci mass appears as 

H 2 3 

M R ~^ (5.4) 

where the gravitational constant G has been introduced, to have correct units. 
Taking again the Sun as the central object and using the previous values of 
Ms, H and i?s, one obtains Mr s» 2.8 Kg. We see that by far the basic 
contribution to M(tj, r) comes from the Weyl mass. However, the situation 
changes drastically near r = 2m or for r\ — > 0, when the 2nd term from (5.3) 
dominates. It is clear that "to be near the singularity r = 2m" or "to wait for 
a long time" have similar effects. 



6 Conclusions 

we studied in this paper a BH (or a compact object with spherical symmetry) 
immersed in a de Sitter space written in a conformally flat form. To be a 
solution of Einstein's equations, we need an imperfect fluid stress tensor on its 
RHS. The proposed model leads to the presence of an inward energy flux that 
diverges at r — 2m. The kinematical quantities associated to a congruence of 
non-geodesic observers are calculated. One finds that the shear and vorticity 
tensors are vanishing but the scalar expansion does not depend on time and 
is asymptotically constant. The cosmic fluid stress tensor does not observe 
the strong energy condition (as for dark energy) due to the negative pressure 
p = —p. We stress the very simple form of the formula (4.10) for the energy flux 
(q asymptotically equals the star energy times the Hubble constant divided by 
its surface area). In addition, it is independent of the gravitational constant. 
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